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DESIGN CHARTS OF STATIC AND ROTARY STABILITY DERIVATIVES 


FOR CROPPED DOUBLE-DELTA WINGS IN 
SUBSONIC COMPRESSIBLE FLOW 

By John E. Lamar 
Langley Research Center 

SUMMARY 

An evaluation of a modified version of the Multhopp subsonic lifting- surface theory 
was made by comparing the theoretical values with experimental data. Near zero lift, the 
theory was found to predict reasonably adequately the lift-curve slope, aerodynamic cen- 
ter, damping in roll, damping in pitch, and lift coefficient due to pitch rate for delta and 
cropped delta planforms and also the lift-curve slope for double -delta planforms. Based 
on this theory, a series of design charts has been prepared for cropped double-delta plan- 
forms in subsonic compressible flow. 

INTRODUCTION 

The National Aeronautics and Space Administration has programs underway to pro- 
vide aerodynamic design information on aircraft configurations and components for speeds 
ranging from low subsonic to hypersonic. Information has already been published that 
indicates the consideration which has been given to both fixed and variable geometry wings 
for use in the design of supersonic transport and military aircraft. (See, for example, 
refs. 1 and 2.) 

A class of fixed wings, the cropped double-delta wings, has found recent application 
in the design of these aircraft. Examples of proposed and actual aircraft using the con- 
cept of cropped double-delta wings are the latest Boeing supersonic transport (ref. 3), 
the Swedish SAAB 35 Draken and SAAB 37 Viggen (ref. 4), and the Lockheed A-ll, also 
designated YF-12A and SR-71 (ref. 4). A search of the literature has indicated the exis- 
tence of some design data, but only a few systematic investigations have been performed 
on cropped double-delta planforms in the subsonic and supersonic speed regimes. (See, 
for example, refs. 5 and 6.) 

In order that a part of this void might be filled, a systematic investigation using the 
modified Multhopp subsonic compressible lifting- surface approach of reference 7 (after 
it was shown to be applicable) with the appropriate boundary conditions was undertaken 



to determine the lift-curve slope, aerodynamic center, damping in roll, damping in pitch, 
and lift coefficient due to pitch rate for nine families of cropped double-delta planforms. 
The purpose of the present paper is to present the results of this investigation in design- 
chart form. The design charts presented are for the attached-flow condition only and do 
not include the effects of leading-edge separation, which are discussed in reference 8. 

SYMBOLS 


A 

b 

Cl 

cl, 


>a 


Cl 


'h> 


'in 


'm r 


c r 


aspect ratio, 


b2 


wing span, feet (meters) 
Lift 


lift coefficient, 


q S 


lift- curve slope, per degree 


Ct lift coefficient due to pitch rate, per radian 


qc 
° 2V 

rolling-moment coefficient, Rolling moment 

’ q Sb 

^OO 

damping-in-roll parameter, , per radian 

a pb’ 

2V 

pitching-moment coefficient about c/4 point, Pitching moment 

q Sc 

damping-in-pitch parameter, per radian 

8 2 £_ 

2V 

root chord 



— _ Mean geometric chord of the total wing 
“ ' b72~ 



Moo subsonic free -stream Mach number 
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m number of span stations where pressure modes are defined 

N number of chordal control points at each of m span stations 

p roll rate, radians/second 

q pitch rate about c/4, radians/second 

q^ free-stream dynamic pressure, pounds/foot 2 (newtons/meter 2 ) 

S total wing area, feet2 (meters 2 ) 

V free-stream velocity, feet/second (meters/second) 

x,y rectangular Cartesian coordinates nondimensionalized with respect to b/2, 

where origin is in plane of symmetry at half root chord (positive x, aft; 
positive y, along right wing panel) 

x a c aerodynamic center, in fractions of c ', measured from leading edge of c' 

- oc 

(positive aft), 21 + - 

’ ac L 4 

y b spanwise location of leading- edge break, feet (meters) 


a angle of attack, degrees 



A outboard leading- edge sweep angle, degrees 

A* = tan - ■'■(tan A//3), degrees 

X overall taper ratio, ?lP- c A°r d . 

Root chord 

X inboard leading-edge sweep angle, degrees 

X' = tan -1 (tan x//3), degrees 
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METHOD OF ANALYSIS 


The method used in this paper for predicting lifting pressures, hereafter called the 
present method, employs the modified Multhopp approach. (See ref. 7.) Briefly, this 
method employs an acceleration potential (developed from a sheet of pressure doublets) 
in conjunction with the linearized Euler equations to relate the pressure difference across 
the wing to the downward velocity over the wing surface. The effects of compressibility 
are accounted for by using the Prandtl-Glauert rule. 

Formulation of the problem of determining the pressure difference across the wing 
leads to an integral equation to which the solution in closed form is difficult, except for 
certain classes of wings, because the answers sought (surface loadings) are a part of the 
integrand. Also adding to the difficulty of solution is the presence of a second- order sin- 
gularity which is in the integrand. For these reasons it has been necessary to resort to 
an approximate solution which makes use of a finite, rather than unlimited, number of 
boundary points over the wing at which the flow is constrained to be tangent. Coupled 
with this finite number of boundary points is a series representation of the lifting pres- 
sures which have the same number of unknowns as there are boundary points and whose 
values are determined by a process of matrix algebra. 

In reference 7 the accuracy of the numerical results of the modified Multhopp 
method was shown to be dependent on the combination of the number of spanwise stations 
m where each of N chordwise control (boundary) points was located. Two methods of 
estimating the appropriate combinations of N and m were studied in reference 7; one 
was concerned with computing the correct leading-edge thrust from the chord loadings 
and the other with seeking the aerodynamic center for which convergence had occurred. 

The converged- aerodynamic- center method was selected to determine the appropri- 
ate pair of N and m to be used in this report because it involved less computational 
labor. Upon examining the results of reference 7 and making use of the method described 
therein, a pattern of N = 8 and m = 23 was postulated to be sufficient for the kinds of 
planforms to be studied herein. In order to verify that this pattern would give results 
which lay in a converged-aerodynamic-center region, the aerodynamic center was calcu- 
lated for several different sets of N and m and is plotted as a function of m in fig- 
ure 1 for one of the families of cropped double-delta planforms for which design charts 
were prepared. A ±2-percent c ' error band is plotted about the N = 8 and m = 23 
results, and for m > 11 the aerodynamic- center values appear to lie within this band. 
However, the predicted results obtained with the N = 8 and m = 23 pattern are not 
necessarily the same as would be found for an infinite number of points, but the two 
answers should compare closely. (The N = 8 and m = 23 control- point pair was also 
used in predicting the rotary aerodynamic characteristics.) 
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Although the present method was developed for steady, longitudinal, lifting prob- 
lems, it can be used to compute rotary stability derivatives because of the relationship 
developed in the lifting- surface theory between tangential-flow boundary conditions and the 
downwash produced by a sheet of pressure doublets. The manner in which computation 
for rotary stability derivatives is made is based on recognizing, as pointed out in refer- 
ence 9, that the steady-state rolling and pitching maneuvers are equivalent to linear twist 
and parabolic camber, respectively. This equivalence can be seen by realizing that the 
upwash on the upgoing wing in steady, rolling flight is just py and the tangent of the angle 
that the section makes with the free stream is just py/V. (The tangent of an angle is 
approximately equal to the angle in radians for small angles.) Integrating this expression 
along the chord gives a straight mean camber line, but one that is twisted along the span. 

A similar procedure can be used to obtain the mean camber line for the pitching motion 
about the reference point. The upwash is just qx for the upgoing part of the wing. The 
tangent of the angle that it makes (small angles of attack) with the free stream is just 
qx/V. Integrating this along the chord results in a parabolic camber. 

VERIFICATION 

In order to verify the applicability of the present method for predicting the static 
and rotary aerodynamic characteristics of the cropped double-delta planforms, a com- 
parison of the results from this method with experimental data (refs. 10 to 18) and with 
results from other theories (refs. 19 and 20) is presented in table I for delta, cropped 
delta, and double-delta wings. It should be noted that in this table no comparisons are 
presented for cropped double-delta wings because experimental data were lacking. How- 
ever, comparisons of data for the other wings with results from the present method indi- 
cate that this method was applicable to wings with breaks in the leading edge and cropped 
tips and, hence, to cropped double-delta wings. Results from the present method are also 
compared with experimental data in figure 2. The theoretical methods selected for com- 
parison are those which predict only the rotary stability derivatives since the theoretical 
methods which predict the static stability derivatives were compared with results from 
the present method in reference 7. In examining table I, it should be remembered that 
the theoretical method of reference 19 has an inherent restriction associated with it 
because its derivation was based on wings with vanishing aspect ratios . 

Static Stability Derivatives 

Generally good agreement is found between the theoretical and experimental values 
for Cl^, where the theoretical values have variations of about ±10 percent from the 
experimental values. (See fig. 2.) Also, the predicted x ac values generally agree 
well with the experimental values and have variations of about ±10 percent from the 
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experimental values or ±5 percent of c. No comparison between the theoretical and 
experimental location of x a c could be made for the double-delta wings of A =1.60 
because reference 18 reports that the pitching-moment data were questionable, and hence 
not presented, for small angles of attack. 


Rotary Stability Derivatives 


Comparison of results from the present method with experimental data and results 
from the other theories for the rotary stability derivatives is restricted to the delta and 
cropped delta planforms since (1) no experimental rotary derivatives were found for the 
double-delta wings and (2) the theoretical methods of references 19 and 20 were not read- 
ily (if at all) applicable to the double-delta wings. 

Damping in roll (C^.- The present method generally predicts values of 

about 20 percent more negative than those measured experimentally for a wing of 
A ^ 1.0. A large part of this loss in theoretical damping may be due to the wing tips of 
the experimental model experiencing a stall condition for which the present method does 
not account. The actual individual differences between the experimental data and values 
from the present method are smaller for all configurations than those resulting from 
use of the method of reference 19, which was developed for slender wings only. How- 
ever, when experimental data and values from the present method are compared with 
values from the method of reference 20, it is found that the predicted values of refer- 
ence 20 agree almost as well with experimental data as values from the present method 
for the delta wings and agree slightly better with experimental data for the cropped delta 
wings. Hence, either the method of reference 20 or the present method could be applied 
with almost equal accuracy for delta or cropped delta wings but the present method is 
more easily applied to cropped double -delta wings. 


Damping in pitch (C m qj.- The present method generally predicts values of C nl q 

which are about 30 percent less stable than the experimental values; however, the results 
obtained by this method agree better with experimental data than those obtained by the 
method of reference 19, which predicts results that are consistently more stable than 
experimental data and results from the present method. The method of reference 19 is 
based on slender -wing theory and is obviously outside of its range of applicability. 


Lift coefficient due to pitch rate (Cl^V- The present method generally predicts 


values of 



about 25 percent larger than the experimental values; however, the 


results from this method agree better with experimental data than the results predicted 
from the method of reference 19. The method of reference 19 is again obviously outside 
its range of applicability. 
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General observation .- If a consistent set of experimental data and the corresponding 
predicted values from the present method, as given in table I, were plotted against aspect 
ratio, it would be found that the correct trends are predicted by the present method for 
all three rotary stability derivatives. 

DATA PRESENTATION 

Because of the multiplicity of planforms that can be encountered in connection with 
composite wings, such as the cropped double-delta planform studied here, design charts 
cannot be expected to eliminate the need for computer calculations to determine accurately 
the aerodynamic characteristics of any particular planform. However, design charts 
which broadly cover the range of planform parameters of general interest for aircraft 
which must fly at both subsonic and supersonic speeds are felt to be useful in selecting 
the general type of planform which might best satisfy the various aerodynamic require- 
ments for the preliminary design thereby reducing the computer study requirements. 
Design charts, which illustrate the degree to which some of the more important aerody- 
namic parameters are affected by the various geometric variables, are presented herein. 

In considering composite planforms, a variety of geometry controlling parameters 
are available for selection as the independent variables. However, since all of the con- 
figurations studied have straight unswept trailing edges, the following variables are the 
ones selected for use herein: 

(1) Inboard leading-edge sweep angle (y) 

(2) Outboard leading-edge sweep angle (A) 

(3) Spanwise location of the leading- edge break y^ 

(4) Overall taper ratio (x) 

These variables, and in particular y and A, were used to develop design charts at 
Moo = 0. However, the charts can be used for any subsonic Mach number if the wing 
geometric properties are corrected in accordance with the Prandtl-Glauert rule. Hence, 
the design charts are presented in terms of these corrected variables y' and A' and 
also in terms of y^ and X. 

Wing geometric characteristics .- Unfortunately, the design charts based on the four 
variables just discussed do not allow convenient identification of other dependent parame- 
ters of interest, such as aspect ratio and mean geometric chord. Therefore, charts of 
wing geometric characteristics have been prepared which show the variation of inboard 
sweep angles and mean geometric chord with aspect ratio; these are presented as fig- 
ures 3 and 4, respectively. Since the wing trailing edge is straight, a knowledge of the 
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mean geometric chord allows a ready determination of the moment reference point c '/4 
with respect to the trailing or leading edge of the root chord. 

Wing aerodynamic characteristics.- The static aerodynamic characteristics for 
the cropped double-delta planforms are presented in figures 5 and 6. Figure 5 shows 
the lift-curve slope and figure 6 the aerodynamic center. The rotary aerodynamic char- 
acteristics for the cropped double-delta planforms are shown in figures 7, 8, and 9. 
Damping in roll and damping in pitch are shown in figures 7 and 8, respectively, and 
lift coefficient due to pitch rate is presented in figure 9. 

Figure 10 is a summary of data which presents the aerodynamic -center locations 
for cropped delta wings as a function of X. This figure is included herein because it 
shows an interesting occurrence; that is, regardless of the leading-edge sweep angle, a 
cropped delta wing has its aerodynamic center at approximately 25 percent of the plan- 
form mean geometric chord when the taper ratio is 0.32. Reference 21 reports a sim- 
ilar finding except that the taper ratio required was slightly different. 

The reasoning used in reference 21 was that if planforms with limiting aspect ratios 
of zero and infinity both had their aerodynamic centers at the 25-percent mean- geometric- 
chord location, then perhaps, by suitable planform variation, the aerodynamic center could 
be kept at the 25-percent mean- geometric- chord point within these two aspect ratio limits. 
By using this basic reasoning, which does not hold for a delta wing (see ref. 22), the taper 
ratio required for a cropped delta wing was determined to be -0.303. 

CONCLUDING REMARKS 

An evaluation of a modified version of the Multhopp subsonic lifting- surface theory 
was made by comparing the theoretical values with experimental data. Near zero lift, 
the theory was found to predict reasonably adequately the lift-curve slope, aerodynamic 
center, damping in roll, damping in pitch, and lift coefficient due to pitch rate for delta 
and cropped delta planforms and also the lift- curve slope for double-delta planforms. 
Based on this theory, a series of design charts has been prepared for cropped double- 
delta planforms in subsonic compressible flow. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, November 21, 1969. 
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TABLE L- WINGS AND REFERENCES USED IN THE VERIFYING PROCEDURE AT LOW SUBSONIC MACH NUMBERS 



Geometric parameters 


Static stability 
derivatives 

Rotary stability 
derivatives 

References 

A 

x. 

deg 

A, 

deg 

572 

X 

Ct 

x ac 

% 

% 

Cl q 







Delta wing 




0.25 

_ 

86.4 



0.006 

0.466 




Unpublished experiment 












.0065 

.459 


— 

— 

Present method 

.50 

0 

82.9 

0 

0 

.0125 

.464 

-0.01 


— 

10 






.0125 

.443 

-.046 

-.049 

-.035 

— 

— 

Present method 
Method of 19 
10 

1.00 

0 

76 

0 

0 

.024 

.022 

.400 

.402 


— 
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.023 

.421 

15 














.023 

.414 

-.087 

— 

— 

Present method 











Method of 19 











1.07 

0 

75 

0 

0 

.024 

.450 

-.075 

-0.72 

1.39 

12, 14 






.024 

.410 

-.093 

-.720 

1.766 

Present method 







— 

-.105 

-1.261 

2.521 

Method of 19 

1.67 


67.4 



.034 

.391 




15 









.034 

.386 

Present method 









2.00 

0 

63.4 

0 

0 

.0425 

.037 

.346 

.375 

-.117 


— - 

10 

11 














.0385 

.375 

-.153 

-.196 


— 

Present method 
Method of 19 










-.161 

Method of 20 










2.31 

0 

60 

0 

0 

.042 

.400 

-.17 

-1.00 

2.05 

12, 14, 16 






.042 

.368 

-.171 

-.995 

2.842 

Present method 







— 

-.227 

-.180 

-2.721 

5.443 

Method of 19 
Method of 20 










2.5 


58 



.045 

.374 




15 









.045 

.363 

Present method 









3.00 

0 

53 

0 

0 

.050 

.048 

.320 

.355 

-.168 



10 

11 














.050 

.352 

-.204 

-.295 


.... 

Present method 
Method of 19 










-.221 

Method of 20 










4.00 

0 

45 

0 

0 

.057 

.362 

-.22 

-1.50 

2.85 

12, 14, 16 






.058 

.338 

-.241 

-1.106 

3.634 

Present method 







— 

-.393 

-.278 

-4.712 

9.425 

Method of 19 
Method of 20 















Cropped delta wings 




0.54 

o 

75.9 

0 

0.303 

0.016 

0.289 




15 









.014 

.256 

Present method 









.89 

o 

67.4 

0 

.303 

.025 

.282 




15 









.023 

.257 

Present method 









.99 

0 

45 

0 

.58 

.0175 

.220 

-0.085 

-0.800 

1.18 

12, 14 






.0255 

.174 

-.098 

-.095 

-.491 

1.727 

Present method 
Method of 20 










1.33 

o 

66.8 

0 

.125 

.030 

.031 

.367 

.338 




13 

Present method 















1.33 

0 

61 

0 

.250 

.032 

.322 




13 









.032 

.279 

Present method 









1.33 

o 

45 

0 

.500 

.033 

.032 

.267 

.203 




13 

Present method 














1.34 

o 

58 

0 

.303 

.034 

.264 




15 









.032 

.258 

Present method 









1.99 

0 

45 

0 

.36 

.036 

.250 

-.16 

-1.09 

2.15 

12, 14 






.043 

.246 

-.185 

-.179 

-.708 

2.709 

Present method 
Method of 20 










3.00 

0 

45 

0 

.15 

.050 

.340 

-.22 

-1.25 

2.60 

12, 14 






.053 

.305 

-.241 

-.239 

-.949 

3.352 

Present method 
Method of 20 















Double-delta wing 




1.30 

82 

60 

0.500 

0 

0.024 

.025 

0.647 

.583 




17 

Present method 













1.60 

75 

65 

.474 

0 

.0365 

.032 





18 

Present method 


.455 












1.60 

80 

65 

.320 

0 

.036 

.032 





18 

Present method 



.488 
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2% c'band 
about the N=Q 
m - 23 point 








(a) Different transformed leading-edge inboard sweep angles are used. A 1 = 50° : A = 0.10; = 0.50. 

Figure 1.- Effect of number of chordwise and spanwise stations on the aerodynamic center for typical cropped double-delta planforms. 


12 




* 2% c' bond 
about the N=8 
m = 23 point 





(b) Different leading-edge break ratios are used, x' = 85°; A 1 = 50°; A = 0.10. 
Figure 1.- Concluded. 




(a) Static derivatives. 

Figure 2.- Prediction accuracy of modified Multhopp method at N = 8 and m = 23 for some aerodynamic characteristics 

of delta and delta-like wings. 
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(c) A' = 50°; A = 0.50. 
Figure 6.- Continued. 
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<h) A' = 72°; X = 0.25. 
Figure 6.- Continued. 
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(e) A 1 = 60°; A = 0.25. 
Figure 9.- Continued. 
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(i) A' = 72°: A = 0.50. 


Figure 9.- Concluded. 
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Figure 10.- Aerodynamic center for cropped delta planforms. 
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